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ABSTRACT 


Simultaneous measurements of waves, orthogonal water 
particle velocities, eeeaturc., conductivity and acoustic 
phase and amplitude measurements at 20 and 40 kHz over 
horizontal and vertical, 2m paths were made at the NUC 
oceanographic tower, San Diego in April 1974. 

Four stochastic models for acoustic propagation including 
geometric, Born, Rytov and Debye approximation models were 
tested. The Krosil'nikov - Obhukov formulation of the Born 
approximation ylelded reasonable results. Difficulty in 
establishing scales obviated good results from the Rytov 
and Debye models. 

Taylor thermal microscales from spatial correlation 
analyses determined vertically across the thermocline varied 
from 6.8 to 178 cm. The measured variance of the index of 
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refraction ranged from 4.5x10- 
Temporal stochastic analyses showed horizontal acoustic 

path losses were ten times vertical losses because of 

diffraction effects. The primary mechanisms for acoustic 


variability were internal and surface waves pumping the 


thermocline up and down. 
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I, INTRODUCTION 


A. HISTORY 

na cONntInuineg pregram of inwsitu ocean microstructure 
research has been conducted at the United States Naval 
Postgraduate School (NPS) since 1971. Under this program, 
a series of three annual experiments dedicated to the study 
of the complex temporal and spatial relationships between 
multiple parameters measured in shallow water have been 
conducted by the Department of Oceanography. Each of these 
experiments has been restricted in Scope to measurements 
inside a 9 m2 volume of water and to oceanographic phenomena 


Whiten Occur in the frequency range 0.001 to 2.5 Hz. 


Bee ObJECTIVES 

Mimeiissthnesis, initial analysis of certain parameters 
measured in the latest of these experiments is reported. 
Primary emphasis is given to stochastic temperature relation- 
ships affecting high frequency (20 kHz) CW sound transmission 
over a short path. The particular objectives of this 
experiment were to investigate ane dependence of acoustic 
phase and amplitude fluctuations on the magnitude and 
eharacteristics of the ocean thermal structure and the 


intensity of the related turbulence. 
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1 ORY 


The description of oceanic turbulence by stochastic 
techniques has followed developments in the atmosphere. 
One of the principal difficulties in application of the 
theoretical methods provided from atmospheric research has 
‘been determination of the range of applicability of certain 
turbulence models in terms of ocean parameters. Another 
- has been the application of these models to sound propagation 
in a manner similar to the relationships developed for 
electromagnetic wave propagation in the atmosphere as 
summarized by Strohbehn [1968]. Historically, attempts to 
resolve these difficulties have been interconnected since 
their primary relevance is to acoustic prediction which 


motivated the research. 


A. STOCHASTIC MODELS FOR THE MEDIUM 

Three expected value techniques which have been applied 
to turbulence calculations in the ocean are 1) the spatial 
correlation (autocorrelation) function, 2) the structure 
function and 3) a "piecewise correlation" function. 
Application of each technique requires a statistically 
homogeneous medium but the order above indicates their 
relative applicability to increasingly less stringent 
Stationarity requirements. The function form for each is 


listed below. 
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The autocovariance function, Rp) : 

Re) = E(t(x)?(xto)] « 1} 
‘The structure function, D,(p) 

D.(p) = E({[T(x) - T(xtp)}°} . (2) 


The piecewise correlation function, RY Cp) 


E{(T(x) - T(xtp') }*} 
Ri Cp) = E[l - SL Cy ee 
a se a 
= aL oa . 2 
o R,(0) 


E indicates expected value and p is the spatial lag. 
POr completely stationary processes, the structure function 


and the "piecewise correlation" function reduce to 


D,(p) 2[R, (0) - R,(e) J ; (4) 


and Ri (Cp) ELL = De euee7 a nenaes) | = fee) /Reae0)' 


(5) 


The stationarity difficulties often encountered in the 
use of Rp) are partly overcome by the use of D,(p) because 
the large scale effects in the ocean are differentially 


filtered out [Black 1965]. The piecewise correlation function 
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can be applied to measurements which have only "piecewise 
stationarity." This makes it the most powerful tool of the 
three techniques. "Piecewise stationarity" is defined as 
stationarity (in the ordinary sense) over a sufficient 
record length for which the structure function and variance 
can be calculated. This method was proposed by Whitmarsh, 
Skudrzyk and Urick [1957}. 

Various analytical equations have been used to approxi- 
mate the form of the auto covariance function used to describe 
the environment. Liebermann [1951] used a simple exponential 
form in his early experiments. Others have used a Gaussian 
form with some success [Stone and Mintzer 1953]. Use of the 
exponential form yields the best fit to experimental data 
above certain spatial lags, but it underpredicts the number 
of "small" turbulent blobs. This undesirable effect is 
eliminated by the use of the Gaussian form. 

Kolmogorov [1941] predicted a 2/3 power law dependence 
LOtemmMercumuegure function for turbulent flow 2answoich there 
was a decoupling of the energy producing and energy dissipa- 
tion regions. Batchelor [1953] transformed this result and 
demonstrated that in wave number space energy spectra have 
a wave number to the -5/3 power law dependence. This was 
first verified convincingly by Grant et al. [1962] for ocean 
measurements and Pond et al. [1963] for atmospheric measure- 
ments. Ocean measurements do not always show the dependence 
since the ocean tends to be stratified much more intricately 


structured than the atmosphere which makes the recording of 


IG: 








continuous stationary data unlikely. However the -5/3 

power law depends not so much upon isotropy as upon the 

clear separation of scales of energy production and energy 
dissipation in which a cascade of energy from larger to 
smaller Seales occurs. Whether or @met suetea form is 
ubiquitous in the ocean is not as significant as the presence 
of horizontal sheets which could cause multiple reflection 
layers. Very sharp, multiple thermal layers have been 
reported in the Arctic [Denner 1971], for example. 

Recently it has been shown that surface wave effects can 
be removed from ocean turbulenee measurements with the result 
that purely turbulent flow remains and a -—5/3 power relation- 
ship may be observed [Thornton and Boston 1974]. These 
findings were based upon measurements both in and around the 
thermocline and showed that anisotropy was not a critical 
difficulty, provided additional energy Sinks or sources did 
not result (e.g. in the thermocline). 

Unfortunately, other energy contributions to the 
turbulence spectra do exist. Microstructure related 
turbulence from salt-fingering is one form [Gregg 1973, 
Denner 1971]. Billow turbulence resulting from shear 
instabilities related to internal waves has also been 
reported by a number of investigators [Woeds and Wiley 1972, 
Gregg 1973, Powell 1974]. Further, the existence of the 
internal wave itself could be expected to #dd energy at low 
frequencies as could surface wave interaction with the 


thermocline. 


14 





B. WAVE PROPAGATION MODELS 

Four methods derived from the Helmholtz equation have 
been used to describe wave propagation through the turbulent 
ocean. These are: 1) geometrical optics, 2) the Born 
approximation [Born 1933], 3) Rytov's method [Rytov 1937], 
and 4) the Bebye approximation [Frisch 1968]. While it is 
not a necessary assumption, most applications of these 
methods presume that the ocean is homogeneous so that the 
scalar wave equation can be used. Such an assumption places 
a rather severe restriction on the acoustic wavelength, i, 


such that 
h << &. (6) 


where = inner scale of the turbulence. This restriction 
will be discussed in detail later. 
1. Geometrical Optics 
The geometrical optics (ray theory) approach is a 
linearization of the wave equation which assumes only 
Rayleigh scattering and refraction effects are important. 
Since diffraction effects are ignored, the high frequency 


(short wavelength) restriction, 
2 
nS TAR (7) 


where L = acoustic path length, is placed on the results. 


Is 








One of the first experiments which used a stochastic 
technique to examine the acoustical effects of small-scale 
temperature inhomogeneities in the ocean was conducted by 
Liebermann [1951]. Through the use of the classical optical 
model, Liebermann was able to predict certain scintillaticns 


in sound amplitude observed at sea by the equation: 


2 _ 8 2 3 
V+ = TS YT © (L/a) (8) 
Oo AI,2 
O 
where ve = variance of the sound intensity, Ty = mean sound 


intensity, Al = the instantaneous deviation of sound intensity 
frem the mean, a? = variance of the index of refraction, 
L = acoustic path length and a = Taylor temperature micro- 
scale [Lumley and Panofsky 1964]. Liebermann also showed 
that sound reverberation encountered in the ocean was larger. 
than that predicted due to thermal inhomogeneities alone. 
This was substantiated by later developments which have 
shown that biological scattering, chemical relaxation and 
bubbles in the near-surface region are important contributors 
ToOnmscound avLtenuan lon. 

Sagar [1955,1957,1959] was able to show that sonar 
amplitude fluctuations resulting from temperature measure- 


ments could be correlated to sea state. Additionally, he 
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showed that the linear dependence of the variance of sound 
pressure level on range [Sheehy 1950] was apparent only when 
thermal microstructure was present. 

A modified geometrical optics technique has recently 
been proposed by Tatarski [1967]. The method is similar to 
the classical one except the path length is described as 
the eikonal or total path length, 98. However, one would 
assume that this would modify one restriction on the results 


Co 
5 | 
\ << L. /0 (11) 


No application of this method to the ocean has yet been made. 
2. Born Approximation 

The Born approximation is a perturbation technique 
applied to the wave equation which is valid only as long 
as the acoustic amplitude fluctuations remain small. This 
is not a serious restriction in the ocean, even under the 
most severe scattering conditions, but the model only accounts 
for a single scattering from a turbulent blob over the 


eouarle, Dauh length, L. This requires that 
A << a vaL (12) 


The geometrical optics results are a special, high frequency 


case included in this model. 


ey 





Stone fia Mintzer [1962, 1965] conducted a series 
of laboratory experiments from which they demonstrated 
solutions to the Born Model which predict the variance of 
normalized sound pressure level, ial in a low frequency 
region (wave region) and in the ray theory region. For 


the wave region they Showed: 


ve = 5 WT k a a L (Gaussian R,(p) (13) 

We = 6 k.* a L (Exponential R,(p ) fale) 

Wo = ELGG )*1 (15) 
O 

X << 2ta (16) 

X << 2nrLh (a) 

Ky = On/r (18) 


Where Es = mean sound pressure level and AP = pressure level 
difference from the mean. For the ray region, no solution 

is possible if the exponential form of the spatial correlation 
function is used since the solution is divergent. But they 
showed that the following equation results from use of the 


Gaussian correlation function: 
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vy? = “= Va aha)? (19) 


Neubert and Lumley [1970] used data from the above 
experiments to compare with atmospheric data based upon the 
Saem model shear nesumesesnowed chat the turbulent scales 
for air and water are of the same order of magnitude for 
"Similar" turbulent flows. They were also able to better 


define the limit on wave length specified in relation (6), 


such that 
L 
ds a/5[o,(u/e,) Ie (20) 
o 2 
a Xp /Oy, (21) 
2 4 
R anh 22 
: a (n/a) (22) 
where: 
de = Taylor turbulence microscale, 
O, = Prandtl number for the medium, 
wees lecall turbutenc velocity, 
Cc, = mean speed of sound in the medium, 
Ry = Richardson number, 
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ne = Kolmogorov turbulence microscale. 
In addition, they found a restriction on acoustic frequency 
Was necessary, though it was usually less restrictive than 


eo. This restriction is 


lbgeel 
22D =a 
em oOLe ae) | )d, Geo) (23) 
It follows from (6) and (20) that 
1 
£, = Ofa/50[o,(u/e,)I°} Com) (24) 


which lends additional meaning to relations (6), (7) and 
er). 

Whitmarsh, Skudrzyk and Urick [1957] conducted 
experiments in the ocean which largely agree with Stone's 
anageMantzer’s laboratory results. However, they noted the 
variance in sound pressure level due to thermal SCRE ECHne 
varied directly as the acoustic frequency and as the range 
squared in the Kolmogorov turbulence range. This result is 
obviously in the intermediate range between the ray and 
Wave solution regions suggested by Stone and Mintzer. 

They also observed that the turbulent patches were 
not spherical, but were spheroidal in shape. In this case, 


a simple transform to cylindrical coordinates shows that the 


a0 





Taylor microscale, a, must be replaced such that 
a® = (b“/3)(2 + c°/b*) (25) 


where b = narrower dimension and ec = longer dimension of 
the turbulent patch. This has the effect of reducing the 
focusing length of an individual turbulent "lens" which, © 
in turn, increases the effective ray solution region. 

The following expression for the variance of the 
Semnd pressure level, ie. has been presented as valid over 


the entire wave and ray regions [Krasil'nikov and Obhukov, 


1956[: 


y" = on a” a L{l - 1/Dtan "DJ (26) 

D= TL/k a. (27) 
For this result, the wave region is defined when 

LES Bee 7 | (28) 


and the ray region is defined for 


A << Ase /p,, (29) 


el 





3. Rytov's Method 

Rytov's method (method of smooth perturbations) is 
the method used by Tatarski [1961] to develop a stochastic 
wave model for the atmosphere. It includes the geometrical 
optics and Born approximation methods as special cases and 
imeludes multiple scattering effects as well. Unfortunately, 
the scattering effects are included such that they cannot 
be parameterized directly in terms of the turbulence. 

Either the spatial correlation or structure function 
may be used with this method as dictated by considerations 
regarding the stationarity of measurements to be analyzed. 


The restriction 


hi << Jee °/L (30) 


is necessary for reasonable solutions of the aed Gant 
equations. It is also necessary to compute the range of 
validity of the perturbation expansion when this method is 
used. This calculation will not be discussed, but Strohbehn 
[1968] suggested that the method is valid as long as the 
variance in the normalized log-amplitude of the Sound pressure 
level is less than 0.8. He also argued that the results 
from sound phase calculations should be less restricted than 
those for amplitude. 

Black [1965] applied this model to certain data from 


imemcuniN experiment. One important conclusion he reached 
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was that the standard deviation of the fluctuation in sound 
intensity through a turbulent medtun, Og» and the principal 


period, T,, of its fluctuation can be predicted by the equations: 


: 3 ee 
4 12 12 . 
Cel opigiacammbn =: 1X (oa) 
ae TYAL/2v_, (32) 
2 2 2 2 2 2 2 2 
Con = Cy Ve ve /10008S + Cy /C, = Crp Ze (33) 
2 
“nm = Dq(L) /b> (33a) 
2 
C.° = D,(L)/L° (33b) 
ae 
rau 3 
eee (36) 
where: 
ys absolute value of the medium's mean velocity 
across the acoustic path, 
T, = mean absolute temperature Ce). 
S = mean salinity in parts per thousand, 
C, = mean speed sound (m/sec), 
Dp (L) = value of the temperature structure function 


at range L, 


D(L) = value of the salinity structure function 
at range L, 


D,(L) = value of the turbulence structure function 
at range L. 


The values for C.. Which Black observed were: 


n 
In thermocline: 5.2 x 107 m7 
Near surface: 120 = 107 m7 
Below thermocline: eeahoe.: Oe 
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4. Debye Approximation 
The Debye approximation to the wave equation is the 
most general model currently available for the description 
of the turbulent ocean effects on sound transmission [Neubert 
1970]. It contains Rytov's method as part of its solution 
and it reduces to the Born model when the single scattering 


condition 
A >> 2d aval (34) 
tseappl1ed. The resultant equation from the model, 


Vmeecuce see i U (35) 


where Le = the Eulerian integral scale (e.g. the largest 
turbulent scale), has no special long range nor high frequency 
restrictions. In other words, it is good when the cumulative 
erfects of small, continuous fluctuations of the local 
turbulent velocity strongly affect phase at longer ranges. 


No use has been made of this model in the ocean to date. 


C. MODEL APPLICATION 

Stewart and Grant [1962] conducted a turbulence experiment 
near the surface in the presence of large waves. They showed 
the Taylor microscale, a, was on the order of 6 cm which 
corresponds to the results (6.4 cm) Neubert and Lumley [1970] 
calculated from Sheehy's [1950] data. Since these measure- 
mencs cOrrespond well and fit the conditions of the current 
NPS experiment, the Neunere and Lumley values were used to 
Provide & rough basis for further analysis. These values, 


coupled with experimental parameters, are shown in Table l. 
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The limits of applicability of each of the wave models and 
the associated value of ve Calculated on the basis of this 
data is shown in Table 2. It would appear from the latter 
table that the Debye approximation and perhaps the Born 
moge for the wave region (if the limit imposed by relation 
(6) is not too strong) would be most applicable in this 
case. The large difference in a resultant from the Born 
model compared to that from the Debye model could be the 
result of poor definition of the micro and integral scales. 
One should also note that the acoustic near-field 
criteria for a piston source imposes an additional requirement 


on the wavelength: 
2 
ees Seve (36) 


where r = radius of the transmitter active disc. Since 
relation (36) requires A > 0.6 om, geometrical optics and 


Rytov's method can be immediately discarded. 


eo 





TABLE 1 


PARAMETER VALUES FOR MODEL TESTING 


Parameter Assumed Value 
Co 1505 m/sec 
v/C, 7 x 107° 
Le 4Om 
de 6.4 em 
Ry 6300 
oO) if 
Roa 5.0 x 10°! 
L 200 cm 
fe 20 khz 
r 7.5 cm 
- ~j 
Ch des exe LO im 
u 10.5 cm/sec 
258 
Ne 58 cm 
a 2.4 em 
Lo 2.2 cm 
O HG. x 107" 
k yb eae 
‘ : 
D ILS 
r 10.88 em 


NL 


Neubert and Lumley [1970] 
Sheehy [1950] 

Black [1965] 

Calculated from above data. 
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Experiment 
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Equation 
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TABLE 2 


RESULTS FROM MODEL TESTING 


Model te Most Restrictive Limit 
=? -2 
Geometric 5.8 x10 ~ (19) A << 1.1 x 10 “em (7) 
Born 6.3.x 10° (13) eo: on << eens) 


(Gaussian form) 


Born fGen’ 107? (14) a 0x 107° << << 15em(12,16) 
(Exponential form) 


Born d < 20.9em (18) 
(Neubert Limits) 


Born (K-0) 7.1 x 107° (26) \ >> 2.9 x 107“em (28) 
Rytov 1.1 x 107° (31) XA << 0.5 em (30) 
Debye 2.4 x ioe (35) None 
<= 40m) 

e 

Debye _2 

(es 3 m) iL ic) ve Mem (25), None 


HOMOGENEOUS ASSUMPTION 
A << 2.2cem (6) 
ACOUSTIC FAR-FIELD CRITERION 
4 > 0.6 em (36) 


Note: Numbers in parenthesis refer to applicable equation. 
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IIIT. EXPERIMENT AND DATA ANALYSIS 


A. EXPERIMENT 

Measurements were taken from May 24 to 26, 1974 from 
the Naval Undersea Research and Development Center's (NUC) 
Oceanographic research tower located 1.6 km off Mission 
Beach San Diego, California. The tower is a stable, three- 
tiered concrete and steel structure built into sandy bottom 
in water of 18.5 m mean depth. It provides enclosed spaces 
for electronic equipment and open-deck working areas for 
handling instruments. Winch controlled trolleys are avail- 
able on three sides of the structure on which devices may 
be lowered to within 1m cf the bottom. Instruments for 
this experiment were installed on an aluminum and steel 
frame attached to a trolley on the west side of the tower 
which faced the predominant swell. A schematic view of 
the experiment layout is shown in Figure l. 

Basic meteorological and oceanographic data were avail- 
able from equipment permanently installed on the tower. A 
bathythermograph trace eyeaimable from the north side of 
the structure facilitated placement of the instrument frame 
in the desired vertical position with respect to the 
thermocline. 

The general environmental conditions encountered were 
moderate winds, sea and swell from the West. Light winds 


to 5 kts, overcast skies and swell to 1 m were usual in the 
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morning (see Table 3). In the afternoon, winds increased 
to 15 kts and seas of 0.75 m with accompanying swell of 1m 
was common. The diurnal variation in thermocline depth was 
approximately le m. 

1. Instruments 

Temperature measurements were made with Fenwell 
Electronics Inc. type k496, double-bead thermistors which 
exhibit low self-heating and well-matched, high-resolution 
characteristics. Calibration of the thermistors, together 
with associated circuitry, was conducted in the laboratory 
Unger closely controlled, stable conditions. A Hewlett-— 
Packard HP-2801A quartz thermometer with an accuracy of 
+ .006°C was used as a temperature Standard and a decade 
box was used as a resistance standard. The same decade box 
Was used to calibrate the thermistor array in the field to 
eliminate resolution difficulties resulting from the poor 
overall accuracy of the decade box. 

Water particle velocities were measured using Marsh- 
McBirney electromagnetic flowmeters (models 721 and 722). 
Since these flowmeters were Paaeiowated by the manufacturer 
only for steady flow conditions, each probe was recalibrated 
for unsteady flow by oscillating it in a water tank at rates 
corresponding to expected experiment velocities. 

A Baylor Co. wave staff system, model 13528R, was 


used to meaSure sea surface elevation. Calibration was 
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accomplished at the NUC tower by successively shorting out 
Sections of the sensor active element while measuring output 
voltage. 

| A Naval Research Laboratory, Underwater Sound 
Reference Division (USRD) type F-33 general purpose uni- 
directional transducer was used as a CW sound transmitter. 
It features two piezoelectric ceramic arrays mounted coaxially 
which, when excited in parallel, exhibit exceptional direc- 
tivity. This was desired to reduce interference from surface, 
bottom and equipment frame reflections. The beam widths 
(-3db down) for 20kHz and 40kHz were 20° and 10° respectively. 
USRD type F-50 lead zirconate-titanate omnidirectional 
transducers were used as receivers. 

The sonar array was calibrated in an anechoic test 
tank filled with fresh water and the hydrophone calibrations 
compared satisfactorily with calibrations previously per- 
formed at USRD.. The frequency standard used was a Hewlett- 
Packard HP-2801A quartz thermometer configured as a frequency 
meter and the voltage standards were an RMS calibrated Fluke 
digital voltmeter and a Hewlett-Packard model HP 3400A RMS 
voltmeter. Phase measurements were conducted using a Dranetz 
model 305 phase meter. Calibration and accuracy measurements 
for the peripheral equipment associated with the sound arrays 
were taken concurrently with transducer calibration. 

The instrumentation accuracy, response and limitations 


are summarized in Table 4. 
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2. Instrument Arrangement 


Two arrays of thermistors were mounted on one side 
of the equipment frame in a common vertical plane. Seven 
thermistors were mounted along a horizontal member and fe 
were mounted along a vertical member as shown in Figure 2. 
The Meeronencation and abbreviated notation used in the 
text is given in Table 5. One thermistor in the vertical 
array was common to the horizontal array. The thermistor 
probes were spaced in an approximate geometric progression 
so that a maximum number of spatial combinations could be 
achieved for adequate definition of the spatial correlation 
and structure functions with a limited number of probes. 
The minimum spacing (5.0 em) was determined by the physical 
size of the thermistor mounts. This was compatible with 
the smallest wavelengths which cculd be measured within the 
time constant of the thermistor. Maximum thermistor spacing 
was determined by physical constraints on the length of the 
Support pipe to prevent strumming. 

Three velocimeters were placed so that the horizontal 
and vertical particle velocities could be measured in the 
Same plane as the thermistor arrays. Additionally, they 
were placed in a vertical line directly under the Baylor 
wave gage with the horizontal component measured from the 
west. The temperature and particle velocity devices were 
offset 29 em from the support frame by mounting brackets to 


minimize false turbulence effects generated by the frame. 
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TABLE 5 


CHANNEL ABBREVIATIONS 


Data Channel Abbreviation Urder in Wigagec 
Vertical Temperatures TMV1-TMV7 From Bottom 
Horizontal Temperatures TMH1-TMH7 From West 

Sound Amplitude SAM1-SAM3 From East 
or 
Sound Phase SPH1-SPH2 ldaeteyagh Abc) 
Sea Surface Elevation WAVE - 
Horizontal Flow Velocity FVAY-FVCY From Bottom 
Vertical Flow Velocity FVAX-—FVCX From Bottom 
Conductivity SAL1-SAL2 From West 
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The sonar array was separately mounted on a tripod 
base which pivoted vertically around the main equipment 
frame but which was held rigid by clamps during recording 
periods. This arrangement facilitated horizontal and 
vertical sound path measurement (Figure 2). The centerline 
of the transmitted sound beam was 57 cm away from the 
vertical plane of the thermistor arrays to minimize unwanted 
reileéctions Irom the other instruments. 

Three F-50 receivers were placed along the center 
of the transmitted beam outside the near field at 2.0, 2.3 
and 2.9 m from the transmitter face. The minimum separation 
(30 cm) was selected to maximize the number of wavelengths 
through which the lowest frequency (20kHz) sound would 
travel between receivers and to remain near the lower spatial 
scales of the measured temperature. The maximum sound path 
length of 2.9 m was selected to minimize reflections from 
the frame and mountings. 

3. Pre-recording Signal Processing 

The signals from the wave height and particle velocity 
devices were recorded directly. All other signals were 
processed through a combination of filters, amplifiers and 
de bias removal devices to maximize signal-to-noise (S/N) 
ratio before recording. 

a. Temperature Circuitry 

Bechmotethe NorazZoncal and vertical temperature 


arrayS was connected to a separate temperature array drawer 


oi 





Which contained Wheatstone bridge networks, amplifiers and 
calibration facilities for seven channels (see Figure 3). 
The Wheatstone bridge network provided the capability to 
match the response of each channel. Additionally, balance 
adjustments were provided to adjust the operating point of 
the amplifiers to the mean temperature so that distortion 
during recording was minimized. 
b. wonar Cireuitre 

As shown in the block diagram (Figure 4), signals 
from the F-50 transducers were amplified and then band-pass 
filtered at the operating frequency. The signal was then 
processed through a second-detector-type demodulator. The 
dec bias was removed from the detected signal and the ampli- 
tude modulation (AM) signal was then amplified and low-pass 
ierlTered . 

Dranetz type 305 phase meters were connected 
between the band-pass filter outputs from receiver hydrophones 
ene and Cwermand, similarly between hydrophomes two and three. 
The de bias was removed from the output of each phase meter 
and the signal was then amplified and low-pass filtered. 

Before each run the de bias setting was made and 
the rms voltage from each band-pass filter output and mean 
phase shift read from each phase meter was recorded. This 
technique allowed much higher signal amplification to be 
used prior to recording because the signal amplifiers were 


not saturated by the large dc segment of the signal. 
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FIGURE 3. BASIC TEMPERATURE MEASUREMENT 


CIRCUIT BLOCK DIAGRAM. ONLY | OF THE 7 CHANNELS IS 
SHOWN IN THE TEMPERATURE ARRAY DRAWER FOR CLARITY. 
TWO IDENTICAL ARRAY DRAWERS WERE USED DURING THE 
EXPERIMENT. 
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4, Experimental Procedure 

Twelve data recording runs of 29 to 73 minutes 
duration were taken at various depths above, in and below 
the thermocline. A total of 28 channels of data were 
recorded for each run including 2 channels dedicated to 
experimental conductivity devices [Kane 1974]. Surface 
elevation, vertical and horizontal particle velocities, 
temperature, sound amplitude modulation, sound phase modu- 
lation and conductivity were measured simultaneously. 

Run depths were selected on the basis of a BT trace 
taken just before the run, but the final run depth was 
chosen to optimize the degree of fluctuation in temperatures 
monitored from the vertical thermistor array. The sonar 
array was held in a fixed position during each recording 
run. Runs were made with the sonar array oriented hori- 
zontally, vertically and at an angle with respect to the 
support structure. This allowed data to be obtained at 
various attitudes through and along the thermocline. A 
summary of all twelve runs is given in Table 6. 

A Sabre III and a Sangamo model 3500, both 14 channel 
FM tape recorders, were used to record the analog data. A 
timing pulse was recorded on one channel of each recorder 
at the beginning of each run to provide temporal continuity 
between recorders. 

The analog recordings were digitized using a hybrid 


system including a Scientific Data System model 5000 analog 
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computer and a Xerox Data Systems model 9300 digital computer. 
Digital data preparation and analysis was performed on an 

IBM model 360 computer. The digitizing data rate used was 
25Hz, however 60Hz electronic noise was coupled into the 
signal, so low-pass filters were used to remove it. Addi- 
tionally, noise spikes which existed in some channels were 
removed by a digital editing routine. Data handling capacity 
was improved by reducing the number of samples by five using 
an inverse transform, numerical low-pass filter [Davidson 
1970]. The latter reduced the analysis Nyquist frequency 

to 2.5Hz. The linear trend and remaining dc bias were also 
removed before analysis. 

Spectral analysis was performed on 26.6- minute 
segments of certain data sets. This choice dictated a 
minimum detectable frequency of 0.00625Hz. A standard lag 
time of 5% of the record length was chosen which resulted in 
40 degrees of freedom for each computed spectrum. The 
resultant 80% confidence limits of the chi-squared distri- 
bution were between 0.73 and 1.30 of the measured spectral 
estimates. A Parzen window was applied to the correlation 
Benicelon during the spectra calculations to minimize the 


effect of finite record length on the spectral analysis. 
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IV. EXPERIMENT RESULTS 


During this experiment, the primary factors influencing 
the results were the sea state, the thermocline gradient and 
the position of the instruments with respect to the thermo- 
cline. The sea state remained nearly constant throughout 
the recording periods although the peak of the swell showed 
a gradual shift to higher frequencies. The thermocline 
‘a depth and gradient continuously fluctuated under the influence 
Of internal waves which are characteristic of the waters off 
San Diego during this time of year. In order to resolve 
some of the effects of the thermolcine variability, analyses 
in and above the layer were required. Experiment runs one 
and two were chosen for this purpose. Additionally, run 
12 was examined to check the results since: 1) it included 
data in and above the thermocline in a single record; 2) the 
thermocline gradient was intermediate between the gradient 
encountered in runs one and two; 3) the sonar array was 
angled toward the vertical rather than horizontal during the 
run and 4) the run was conducted two days after runs one 


and two under slightly different sea conditions. 


A. ENVIRONMENT 

The thermocline migrated diurnally from near the surface 
in the morning to about 14m in the evening. It was generally 
fee use near the etn raice but near the bottom it was compressed 


so that a very sharp thermal gradient (2°C/m) was created. 
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Superimposed upon this effect was an apparent four-—-hourly 
dilation and contraction of the isotherms. 

Figures 5-9 show representative bathythermograph (BT) 
traces taken periodically during the experiment. The traces 
are plotted so that the 12°C mark corresponds to the time 
the BT was taken. Recording run envelopes in depth and time 
are included to show the position of the instrument array 
relative to the thermocline. The 11°, 12° and 13°C isotherms 
‘ are shown to indicate the relative intensity of the tempera- 
ture gradient during the runs. Figure 9 is an expanded 
time plot for BT's taken during run 12 which shows a large 
internal wave having a 36 min period; this was observed 
throughout the experiment. 

The four-hourly variation in the thermocline gradient 
probably was caused by internal wave interference effects 
Since analysis showed that these variations were associated 
with periods of high temperature variance but not with sur- 
face waves nor heat budget considerations. 

Surface slicks presumably associated with internal waves 
were seen propagating from the West at certain intervals 
during the experiment. A typical example of the effect of 
these waves is shown in the vertical temperature analog 
records of run number 12 (Fig. 10). The large 36 min wave 
peak passing the vertical temperature array is shown by the 
dramatic temperature decrease which occurred at about Ngaya fe) - 


Peme and 1.5 min period, small-amplitude internal waves 
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effects are seen superimposed upon the larger internal wave. 
The "grass" is surface swell related. Spectral analysis 
showed that a broad band of internal wave periods from 

160 secs (maximum detectable by the method) to as low as 

30 secs were superimposed on the large wave as well. 

Figure 11 shows the resultant sound effects from the 
internal wave temperature perturbations of Fig. 10 during 
run 12. The sound array was angled 25° to the vertical in 
‘run 12. The sound amplitude modulated signals SAM1 and SAM2 
were initially above the layer and were not affected by 
the turbulent thermocline losses as severely as SAM3 (below 
the layer). As the layer moved up, the interference and 
diffraction effects of the turbulent blobs [Whitmarsh and 
others 1957] in the vicinity of SAM3 were diminished. This 
is shown by the increase in sound pressure level (about 
107? wb) on the SAM3 analog record. 

When the layer moved into the vicinity of SAM1 and SAM2, 
the sound pressure level at these transducers was reduced 
by about 6 x 107 ub and 3 x tometer respectively. The 
thickness of the "grass" on these plots is an indication of 
the wave-induced amplitude scintillation effects (~2 to 5 
X 107 | ub) caused by pumping the turbulent blobs back and 
forth across the sound path. (Similar results occurred when 
the sound path was horizontal, but the scintillations were 
then caused by waves pumping the thermocline up and down. ) 
Additionally, the movement of the layer (AT = .3°C) between 
SAM1 and SAM2 caused sound phase modulation, SPH1, to phase 


martc about 30°. 
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B. TEMPORAL STOCHASTIC ANALYSIS 

Time series analyses were performed on 26.67 minute 
Segments of each of the sound channels, the wave record, 
three temperature records and two water particle velocity 
records. Analyses included for each signal the calculation 
of the first and second moments, variance, mean period 
(calculated from total zero upcrossings), autocorrelation 
function and autospectrum. Also, the crosscorrelation, 
quadrature and co-spectra, phase difference, coherence- 
squared and linear regression analysis were calculated 
between specific pairs of channels (see Table 7). Examples 
of cross-Spectral analysis results for runs two and le are 
EHOW ame Pigs. l2-15. Additional results are given in 
Appendix A. Thermistor 2 on the horizontal array, TMH2, 
was selected for temperature comparisons with the sound 
records because this thermistor was most central to the 
sound propagation path during the runs. 

Mi bas Clo babiSlLcs 

Representative basic statistics of temperature, 

waves, sound amplitude and sound phase are shown in Table 
8. By examining the order of magnitude of these statistics 
in relation to the depth of the run and the proximity to 
the thermocline,certain qualitative relationships between 
mae wOaromebers can be seen: 

1) The magnitude of the temperature variance was 


determined primarily by proximity to the layer and its mean 
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TABLE 7 


CHANNELS COMPARED IN TEMPORAL STOCHASTIC ANALYSES 


Abscissas For Regression Analyses 


» Channels TMH2 WAVE SAL2 
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FIGURE 12. Run 2 TMH2 vs WAVE Spectra 
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SAMPLE LENGTH - 26.67 min 
FREQ. RESOLUTION - .00625 Hz 
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period was controlled by surface wave swell until the layer 
was very close to the thermocline. 

2) Sound amplitude variance was controlled by proximity 
to the thermocline (temperature gradient) when the sonar 
array was horizontal. The movement of the thermocline 
with respect to the sound path must therefore be the primary 
influence on sound amplitudes in this configuration. 

3) The amplitude variance when the array was nearly 
vertical was much smaller than when it was horizontal. 

This indicated the effect of turbulence proximity to the 
hydrophone was smaller than the effect of temperature 
gradient motion over the sound path. 

4) Since the mean period of the sound fluctuations was 
much smaller than the mean period of the swell at depth, 
turbulent interactions related to the higher frequency 
internal wave motion of the thermocline were emphasized 
an the sound measurements. 

5) When the sound array was vertically oriented vice 
Merizontal;, much larger phase variations occurred even 
though amplitude variations were smaller. This resulted 
from the large temperature variations over thepath length 
caused by waves pumping the thermal gradient up and down. 
This also showed that the temperature variation effects 
on sound amplitude and phase were in quadrature. 

6) The order of the mean sound phase variances indicates 


that phase was primarily influenced by surface waves when the 
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array was horizontal. Since spectral analysis showed that 

the wave coherence with temperature was very low at the 

surface wave frequencies, this effect may be bubble related. 
2. Spectral Analysis 

Cross-spectral analysis of the temperature, waves 
and sound records for the horizontal sonar runs one and 
two showed that the sound amplitude changes (Figs. 13 and 
14) were in phase with the waves and 180° out of phase with 
‘temperature in the regions of highest coherence. This was 
expected in the stable density structure encountered. 
Additionally, the coherence squared function (simply re- 
ferred to as coherence) showed that only the long-period, 
prominent swell was important to the sound amplitude modu- 
lation because the higher frequency wave particle velocities 
were much more rapidly attenuated with depth. 

Sound phase (Fig. 15) was found to vary directly 
with temperature at the internal wave and surface wave 
frequencies although coherence at the surface wave frequen- 
cies was much smaller than at the internal wave frequencies. 
This indicated that the turbulence associated with the layer 
was more important than bubbles or surface wave perturbaticn 
of the thermocline in influencing sound speed over the 
horizontal sound path. A comparison of the coherences of 
Peco le waned, §S°shews vuhat the temperature effect is at 
least as important as. bubbles even at the surface wave 


frequencies. 
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Similar cross-spectral analyses with the sonar path 
angles 25° from the vertical (run 12) yielded markedly 
different results. For the sound path entirely above the 
layer (SAM1), coherence with temperature waS small over all 
frequencies. Small coherence peaks occurred at 0.15Hz and 
0.7Hz which indicated that the vertical motion of the turbu- 
lent blobs forced by surface swell had some effect in the 
angled sound path. 

Sound phase coherence with temperature was largest 
af the internal wave frequencies but smaller, broad-banded 
coherence at the surface wave frequencies was also observed. 
This confirmed that the phase was being controlled by the 
vertical motion of the thermocline. 

Most of the sound path associated with SAM2 (2.3m) 
Was not in the thermocline during the majority of run le. 
However the hydrophone was near the top of the layer. AS a 
result, its amplitude was in phase with temperature at the 
internal wave frequencies and coherence at other frequencies 
(eco ia veAz COmeeliiiom Veit Lege Ghiau tnesintlertference effect 
of the proximity of the thermocline-associated turbulence 
to the receiver was theprimary effect on SAM2e variability. 

Very high coherence and 180° phase difference between 
TMH2 and SAM3 at the internal wave and the surface wave 
frequencies was observed. Since the acoustic path associated 
with SAM3 generally spanned the thermocline, this indicated 


that the variability of amplitude across the entire 
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thermocline was controlled by refraction and diffraction 
Simigawetoethe horizontal ease.  Themeffjecteot diffraction 
was smaller than for the horizontal case, because the len- 
ticular shape of the turbulent blobs present broader profiles 
to the vertical sound path. The reduced diffraction effects 
caused the variance in sound amplitude for run 12 to be an 
order-of-magnitude less than that for run l. 
3. Regression Analysis 

Sound amplitude and phase variances, v° (SAM) and 
v°(SPH), were plotted versus TMH2 variances, v(t), and 
least squares regression analyses were applied to the data. 
Linear relationships between v°(SAM) and v° (SPH) versus 
v“(T) were found to exist for runs one and two but not for 
run l2. This was because of insufficient variability in 
the portion of run 12 analyzed. The empirical relationship 
below was developed from the mean regressions calculated 
hor w@aums One and two: 


ae (37) 


v°(SAM) = (0.34ub?/°C*)V°(T)e 
where L = acoustic path length and x = the distance required 
to reduce Vv“ (SAM) by 1/e (determined empirically). x is 
near the size of the Taylor microscale. 

The derivation of Equation (37), plots of the 
variances over the run segments and the regression plots 


are given in Appendix B. 
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C. SPATIAL ANALYSIS 

Temperature spatial correlation and structure functions 
were calculated directly from ensemble averages of 26.67 min 
Peeords from all possible combinations of thermistor pairs 
along the vertical thermistor array. Twenty-two spatial 
lags with values between 0 and 260 cm were thus determined 
from seven thermistors. The value of Cr (see Equation 33a) 
and the Taylor microscale as well as the Eulerian integral 
scale were derived from the above functions (Table 9). 
Mowever, since large variations were seen in these functions 
the validity of the calculated scales was questionable. 

In an attempt to remove the non-stationary effects 
resulting from internal waves, analyses using ensemble 
averages of short segments having nearly constant turbulent 
intensity were conducted for runs two and twelve. The 
segments were chosen so that examples of low, medium and high 
turbulent intensity could be examined (Table 10). Because 
the thermocline steepness was highly variable along the 
vertical thermistor array, however, only the medium inten- 
sity correlation function was smooth enough to provide satis- 
factory Eulerian integral scales and the structure function 
of all samples varied widely. Values of the Taylor microscale, 
Fulerian integral scale and an generated from hand-smoothed 
results are given in Table 11. The correlation function and 
Structure function values as well as the mean variance for 


the vertical array thermistors are given in Appendix C. 
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TABLE 9 


Temperature Spatial Scales (cm) for 26.67 min Ensembles 





Run # Taylor? Taylor Eulerian 2 
Microscale Microscale Integral Scale Gp 
mo (15 uo 36 1.45x1074 
2 178 158 175 3.6x107> 
12 60 100 36 8.3x107> 
TABLE 10 


Relative Turbulent Intensity for Special Spatial Analyses 


Qualitative Ensemble TMH2 Mean Range Run 

Intensity Size (min) Variance(°C)¢ (°C)e — 

Low 4.67 0.27x1071 *0.77x10°° «12 

Medium 1.67 4, 2ux1072 +4.95x107+ 2 

High 3.73 7.79x10-7 £2,23x10 7 2 
TABLE 11 


Temperature Spatial Scales (cm) for Special Ensembles 


il 2 


Qualitative Taylor Taylor Eulerian C 2 

Intensity Microscale Microscale Integral Scale ee. 
Low ae 6.8 oy 3.3x107> 
Medium 71 0 35 4 4x1079 
High is 20 63 5.4x107> 


Taylor! = Taylor microscale from spatial correlation 
Taylor® = Taylor microscale from structure function 
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D. APPLICATION OF MODELS 

While it is apparent that some of the scale estimates 
calculated from the vertical thermistor array may be erron- 
eous, it is instructive to apply the theoretical models 
for the sound amplitude variance to see which model is 
most applicable in the non-stationary case. Table l2 lists 
the equations necessary to calculate the basic parameters 
for application of the models. For the purposes of this 
examination, the values of the Richardson number CR.) , the 
Prandtl number (o,) and u/c. were assumed to be the same 
as reported by Neubert [1970] (see Table 13 for values). 
The speed of sound is used to calculated Wilson's Equation 
(Urick 1967]. 

Tables 13 and 14 list the resultant parameters of each 
Sample for which scales were determined. It is noteworthy 
that the values of the refractive index (a) calculated above 
are an order of magnitude above the near surface results 
reported by other investigators [Haley 1972, Medwin 1973] 
(see Table 15) except in the low turbulent intensity cases. 
Miers) is further evidence that the relative influence of the 
muebulent thermocline in sound Scatuering is much greater 
than bubbles. 

If the order of magnitude of the Taylor microscale is 
correct (see Table 13 for comparative values) the values of 


aa are one to four orders of magnitude greater than previously 
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TABLE 12 


EQUATIONS USED FOR MODEL PARAMETER EVALUATION 


Parameter Equation 
Cy Cc. = 1449 + 4.67 — 0.557° 77 0. 7a 
oO am = rh. 6°Ve(T) + 0.055°v"(T) I/e,2 
re de = avo, = 2.65a 


he (R,/15) = 2.07 x 10/a! 


= 
= 
iH 


2 ae 2 
Cs C = Cn Looe 
% % 

L, Ly = 0(2.2/2.4a) = 0(0.9a4) 
Al = mean temperature (°C) 
a = depth (m) 
vo (T) = variance in temperature (°C) 
a = Taylor microscale (cm) 
Ne = Kolmogorov microscale (cm) 


* equation for £. assumes u/c. is a constant 
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ots ll 3 


Experiment Parameter Values for Model Tests 


Parameter RUN 1 RUN 2 RUN 12 

6. (m/sec) 1516 1498 1515 

ao” Ny esol 2.1x107> 7.9x107 

L (cm) 200 200 200 

a (cm) 15 178 60 
a 13.5 160 5h 

A (em) TaD a5 tab 

k, (em™*) 0. 84 0.84 0.8 

Ch 1.45x107" 1.05x107" 2a al ai 

L, (em) 36 175 36 

O1 7 7 7 

Me (cm) OQ 472 159 

Ne (m) OES 21.4 19.2 

R, 6300 6300 6300 

ama 3.0x107° io Qealn 3.7x107? 
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TABLE 14 


Short Sample Parameter Values for Model Tests 


Parameter 


oe (m/sec) 


O 
L (em) 
a (cm) 
x 
fo 
A (em) 
-1 
Ko (cm ~—) 
c 
n 
E (cm) 
°K 
do (cm) 
No (cm) 
Re 
a-a (cm) 


Low 


Intensity 


1512 

EO ahOme 
200 

23 

21 

7.5 
0.84 

1 05x07 


ay 


Medium 


Intensity 


1508 
2.0x107> 
200 

71 

64 

aD 

0.84 


1.16x107* 


High 


Intensity 


1500 
2.7x1073 
200 

15 

13 

feo) 

0.84 
1.29x107 
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reported. This result has a direct relationship (in the 
Born model) to the variance of sound amplitude. 

The variances of sound amplitude calculated by each 
lode) previously discussed (Geometric, Born, Rytov, Debye) 
are presented in Tables 16 - 19. In general, the best fit 
mo. the measured sound variances were calculated from the 
Born model using the Krasil'nikov-Obhukov formulation. A 
Surprisingly good fit to the true variance was also obtained 
from the geometric model in the mid-intensity range. Because 
of the above, the stimates of the Taylor microscale seem 
reasonable. 


Since the accuracy of the scales, C £ and L.> the Eulerian 


E 
integral scale, is in question, no firm statement can be 
made about the adequacy of the Rytov and Debye models. 
However, because L, is thought to be correct for the mid- 
intensity ensemble, the Debye model appears incorrect. 
Presumably, calculation of the Structure function and corre- 
mitlon function from the thermistors in the horizontal array 
over Stationary ensembles would resolve this question. 


Unfortunately, this was not accomplished in time for 


maclusion in this thesis. 


12: 








TABLE 16 


Geometric Model Experimental Results 


ve = 4/15 vu a*(L/a)?> 3 i aou<< ay (cm) 
| 2 2 o 2 
Ensemble Model V Observed V~ (ub) £. /L (em) 
Ron 1 Obese 1.12x107" 0.91 
Run 2 3x107° 1.45x1073 128 
Run 12 1.1x107° 405x107 146 
Low Intensity 7.8x107° 7.1x107° ae 
Mid Intensity 4 .2x107° 9.5x107- . on 
High Intensity 8.2x107> 3.0x107" Oe 3's 
TABEE 17 
Born Model Experimental Results 
[Krasil'nikov and Obhukov ] 
Vo = k,*a°ab(1-1/Dtan>p) ; D-= mL/k a 
Ensemble Model ve Observed v- D 
Run 1 ee Oe 1.12x107" 50 
Run 2 7.9x107" 1.45x107> 4.2 
Run 12 3.74x1072 4.05x1072 125 
Low Intensity 5. Geo” 7. sale’ 32.5 
Mid Intensity » Gear: 9.5x107- 10.5 
High Intensity eel Seo clans 50 


D>> 1 7+ Wave Region 


(es 








TABLE 18 


Rytov Model Experimental Results 


fel 6b 11/6 -7/6 ; 7.5 << 7 


5 
Observed Vv“ wt, /L 


Ensemble Model ve O 
Run 1 5.2x107 > 1.12x107" 5.5 
Run 2 2,.8x107> 1.45x1073 150 
Run 12 iessao = 405x107? 25 
Low Intensity 2.8x1072 7.1x107° 9.9 
Mid Intensity 3.4x107> 9.5x107> by 
High Intensity 4 .2x107> 3.0x107" rae 
TABLE 19 
Debye Model Experimental Results 
Vo = 20°k “LL 
Ensemble Model ye Observed ve 
Run 1 2.1x107> 1.12x107" 
Run 2 0.22 1.45x107> 
Run 12 6.3x10 > 405x107? 
Low Intensity 1.7x107°> 7 isa © 
Mid Intensity 4x107° 9.5x10> 
High Intensity 0.13 3.0x107* 
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V. CONCLUSIONS 


The theories of turbulence which have been developed 
for the ocean require some degree of stationarity in the 
medium before they may be applied. This paper has examined 
four basic stochastic models of acoustic wave propagation 
under the most severe conditions of non-stationarity, e.g. 
near a sharp thermocline which is being oscillated by both 
large surface swell and by internal waves. 

The experimental conditions specified that the solution 
form would be in the wave region for which the Born model 
would be most accurate. The solution of the Krasil'nikov- 
Obhukov form of the Born model from directly measured, 
Oc€an parameters gave reasonable results. 

The solution of the Rytov and Debye models for the same 
conditions yielded results two orders of magnitude higher 
than was actually measured. Because evaluation of scales 
was difficult, this was expected even though both models 
are more general than the Born model and, in fact, contain 
its solution in a subset. This also showed the determination 
of the Taylor microscale was more reliable than the deter- 


fumaclOn Of the structure constant, © 2 or the Eulerian 


At 3 
integral scale for measurements across the thermocline. 
An attempt was made to apply the correlation function 


and structure function piecewise to small sections of data 


Meweaewere thought to be relatively stationary. For 


(e 








application across the thermocline, this method did not work 
Well becauSe Selection of such sections is very difficult. 

The results of the time series analyses showed that Eye 
losses over the acoustic path were greatest when the array 
was horizontal near the thermocline because the diffractive 
losses were an order of magnitude greater than those for 
the nearly vertical case. Interference effects due to the 
proximity of turbulence to the receiving hydrophone accounted 
for a 0.14% signal loss and scintillations due to surface 
Swell accounted for 0.07% signal modulation. The effects 
of bubbles were small because the measurements were rela- 
tively deep and the acoustic frequency was much lower than 
the peak bubble resonant frequency. Measured loss over the 
2m path was on the order of 5% and phase shifts of up to 
30° occurred for these worst-case conditions. 

Current sound propagation models based upon empirical 
data have undoubtedly included these effects in surface 
Scattering or reverberation. But the roughness description 
for the sea surface used by sonar models is based upon sea 
State. The losses here described are, instead, related to 
the swell wave-induced water particle velocities and internal 


waves. 
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APPENDIX A 
Spectral Analysis Plots 


This appendix includes examples of PE OeeA ECO R GTA 
analysis plots so that qualitative comparisons between 
sound, temperature and waves during the experiment can be 
made. Included are spectra for 1) sound amplitude versus 
temperature, 2) sound phase versus temperature, 3) wave- 
related effects and 4) comparison of one sound channel 
with another for the same run. The plots included are 


listed below. 


Sound Amplitude Versus Temperature 


Run # Channels Compared 
als SAM1 - TMH2 
1 SAM2 = TMH2 
1 SAM3 = TMH2 
2 SAM1 = TMH2 (see text) 
2 SAM2 - TMH2 
2 SAM3 = TMH2 
12 SAM1 - TMH2 
le SAM2 = TMH2 
12 SAM3 - TMH2 


sound Phase Versus Temperature 


Run # Channels Compared 
1 SPH1 - TMH2 
2 SPH1 - TMH2 (see text) 
2 SPH2 - TMH2 
2 SPH1 - TMH2 


ai 








Wave-Related Effects 


Run # Channels Compared 
1 TMH2 WAVE 
1 Syge ll WAVE 
2° SAM1 WAVE 
6 TMH2 FVBY 
le TMH2 WAVE 
sound Channel Comparisons 
Run # Channels Compared 
° SAM1 SAM2 
2 SAM1 SAM3 
2 SAM2 SAM3 
2 Sledelal SPH2 
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APPENDIX B 


Analyses of Variance Regressions 


Variances of 56 sec samples of data which included 1400 
data points each were calculated along 26.67 min segments 
of certain channels of runs one, two and 12. Regression 
analyses for sound amplitude and phase versus temperature 
were computed using these variances. The resultant plots 
of the variances with the regression curves are shown after 
Table B5. A summary of the regression equations is given 
in Table Bl. The following mean variance equations were 
computed from the average of the regression equations listed 


for runs One and two: 


V°(SAM1) = 3.2 x 107° V°(?) (B1) 
V°(SAM2) = 1.6 x 1072 Vo(T) (B2) 
v> (SAM3) = Hi x 107" v°(T) (B3) 


ity v°(T) is eliminated from the above equations (stationary 
assumption), then because SAM1 and SAM2 are 30 cm apart and 
60 em separates SAM2 and SAM3: 

ve (SAM) o¢ ($)7"73° om (BY) 


where L = sound path length. 
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This equation can be checked by comparing the sound 
amplitude variances for specific 56 sec samples (where the 
sample is stationary) and similarly deriving the relation- 
ship 


-~L/x (B5) 


V° (SAM) o (5) 
where x is the distance required to reduce v~ (SAM) by Of 
This is done in Tables B3 - B5. The following equations 


result: 


: s-L/34.5 
Run 1: V (SAM)< 5 cm (B6) 
j 
Run 2: V°(SAM)c one em (B7) 


Conversion from base = to base e and placing the path 


length reference at lm vice 2m yields (from equations Bl, 


B4Y, and B5): 
V°(SAM) = (0. 34ub?/°c*)Vo(T)e7L/% (BS) 
where: 
L = acoustic path length in meters, 
x = temperature spatial scale for which 
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V°(SAM) = V°(SAM)_ /e 
v<(T) = temperature variance 
v° (SAM) = variance of sound amplitude. 


And the spatial scales in (B4), (B6é) and (B7) become 43cm, 
50cm and 61lem, respectively. 

Since v° (SAM) is known to be a function of ee the 
variance of the index of refraction, and q = 2.1 x 10°? 


for run two (which had the best regression fits): 
b = 
WS) =e eed kOr das Voile to (B9) 


Assuming that x << L, the following approximation can be 


made: 


Re era (B10) 


x 
yer 9) 5 


Substituting Equation B(10) into Equation B(9) and multi- 
\ 
plying by ie to compensate for the spherical spreading 10Sss 


yields: 


V° (SAM) = 7.7 x 10" ae ve (T) xL (Bia) 


Tou 








which is similar to the exponential form of the Born 


approximation: 


abst Mee Pe otmelartiy Tacluded in Equation (B11): 


mpc 


V°(SAM) amie le =X 107° Ko oO v°(T) Boley: 


The plots following Table B5 are included to allow 


(B12) 


(BT3) 


qualitative comparisons of conditions between runs and of 


the variances between channels during the same run. 


For 


each run analyzed the variance versus time plots are given 


followed by regression plots between various channels. 


In OV 





TABLE Bl. Regression Equations for Sound Versus Temperature 


Variances 
AMPLITUDE 
Channel Run # Equation 
SAMI 1 V°(SAM1) = [3.4V°(T)] x 1073 
2 V°(SAM1) = [3.0V°(T)] x 1072 
1: v°(SAM1) =) (ia ov" (T)]. x 1074 
2 = 2 -3 
SAM2 | 1 V°(SAM2) = [2.0V°(T)] x 10 
2 v°(SAM2) = [1.2V°(T)}] x 1073 
12 V°(SAM2) = [0.7 + 1.5V°(T)] x 107° 
2 _ 2 -i 
SAM3 fl V°(SAM3) = [4.8V°(T)] x 10 
2 V-(SAM3) = [3.3V°(T)] x 107" 
12 v°(SAM3) = [0.2 + 4v2(T)] x 1074 
PHASE 
Channel Run# Equation 
SPH1 1 V°(SPH1) = 1 + 41V°(T) 
2 V°(SPH1) = 0.2 + 0.7V“(T) 
12 V°(SPH1) = 3.7 + 77V°(T) 
SPH2 fl V°(SPH2) = 2 + 8V“(T)# 
2 V° (SPH2) = 2 + 22V7(T) 
12 V°(SPH2) = 2 — 25V“(T) 


*¥From 20 minute record. 
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TABLE Be. 


Run 1 Variances 





samp le Run 1 Run 1-18 Run 1-13 Run 1-25 
ST arice 26 2/3 min 56 sec 56 sec 56 sec 
reer) 2.2) senlan %.1x10° 3.2x10° 0.43 x10° 
V.-. (ub?) TO sO? «S20 on? “Wosoeeao™ 186 gon 

SAML 

ub) 7.6x10°  39.0x10° 7.2x10° 13x10 
Veang(hb*) 2.5 x 10> 9x10°  3.3x10° 3.2x10° 
TABLE B3. Sound Amplitude Distance Loss Factors (Run 1) 

sample Rm 1* Run 1-1 Run 1-2 Run 1-3 

Ricesiice 26 2/3 min 56 sec 56 sec 56 sec 
30 em Ls rhs ee 1.67 1.38 
60 em 3.04 4,33 2 8 4.06 
90 em 4 oO pie ie 3.64 5.63 

quantity mean standard normalized accepted 
loss deviation loss loss 
; factors in loss factors factors 
distance factors 
30 em 146 Oy als 1b - (a0 eee 
60 cm 3.52 eles: Bp Nal 2. 30 
90 em 5 02 1.19 3.44 3,45 


not included in mean loss factors 
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TABLE B4. Run 2 Variances 








Sample Run 2 Run 2-9 Run 2-26 Run 2-20 
26 2/3 min 56 sec 56 sec 56 sec 
variance 
CC) 46.2 x 107° 98.1 x 107° 143.4x 107° 121.4 x 107° 
Vea (ub) 140.x 107° 311x107 422x107 420x107 
eens) 56.0x10-° 122x107? 146x107 178x107 
Bath”) EO RM. 80.5 2 ie Mey ase mee Isa ae? 


TABLE B5. Sound Amplitude Loss Factors with Distance (Run 2) 


Sample Run 2% Run 2-9 Run 2-26 Run 2=20 
Ailecance 26 2/3 min 56 sec 56 sec 56 sec 
30 em Cw ES 2.9 es 
60 em 3.29 4.0 3.62 Cre) 
90 em 8.24 10.2 Gree 7 oval 

quantity mean standard normalized accepted 
loss deviation 1oss loss 
distance factors Tmlocs, factors factors 
factors 
30 em 2.6 Ones i eae 6) 1.4 
60 em B70 O22o Ale) 2.8 
90 ecm 9.9 0.73 350 4.2 





¥not included in mean loss factors. 
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APPENDIX C 


Spatial Correlationwand Structume Punctions 


The values of the spatial correlation function and 
Structure function calculated during spatial stochastic 
analyses are tabulated in Tables Cl - C6. The functions 
were calculated along the vertical thermistor array. 

The variance of the temperature at each thermistor over 
the ensemble is also given. 


The Tables include: 


Table # Ensemble 
CH Run 1, 26.67 min 
C2 Run 2, 26.67 min 
C3 Run 3, 26.67 min 
cy Low turbulence intensity 
C5 Middle turbulence intensity 
C6 High turbulence intensity 
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TABLE C2. Run 2, Vertical Spatial Correlation and 


Structure Function 


Thermistor No. Variance (°0°) 
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Ong 0.3607 5.0000 
1.0090 0.3676 15 00d 
0.9452 0.0160 20.0000 
1.0354 0.0388 25.0000 
1.0891 0 SAGs 30.0000 
0.9706 0.0749 45.0000 
ES 03 OR Shes: Sys) Oe (ON0 
0522 0.2159 80.0000 
1.0844 0.6326 85.0000 
0.9813 0.2662 UNECE CoCo 
0.5463 Oy2oue4 1355.0 000 
0.6164 0.6939 ssl em enenele 
OSS, 0% Sf TD C000 
OS JLe 0.4553 160.0000 
Je Dail 0.9197 I .0000 
0.4985 Oss) Sie 180.0000 
0.5580 0.4413 180.0000 
0.4878 0.6848 205.0000 
0.5458 0.6783 aap. 0000 
O07 Siay 0.9619 260.0000 
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fees C3. Run 22; Vertical Spatial Correlation and 
Structure Function 





Thermistor No. Variance (°C) 
a Ore 
2 0.2433 
3 0.1690 
y 0.1444 
5 0.2275 
6 0.0743 
rf 0.0083 
Correlation S timuet use Lag 
Bumer on Funetion 
1.0000 On 0 G0 
1.0558 020211 5.0000 
a6 4 2 O26 301 15.0000 
1.0859 0.1540 20.0000 
0.7096 0.1640 25.0000 
1.5349 1.0283 25.0000 
0.7463 0.2231 30.0000 
0.7552 0.3984 45.0000 
0.1447 0.3588 55.0000 
0.1852 0.6946 80.0000 
0.1969 0.8039 85.0000 
0.1995 0.9920 100.0000 
0.8057 1. 6228 135.0000 
0.6778 1.6854 E50 0000 
0.6295 1.6494 55.0000 
0.6143 2 ene 160.0000 
0.5148 Bees Bil 175.0000 
0.4773 C.D) 180.0000 
O21 gn 8) 9s 180.0000 
0.2998 2 lee 205.0000 
0.1487 iL SGwo2 235 .0000 
0.1214 2s Oites 260.0000 
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TABLE CH. Low-Intensity, Vertical Spatial Correlation 


and Structure Function 


Thermistor No. Variance (°¢%) 





1 0.0666 
2 0.0774 
5 0.0062 
4 0.0038 
2 0.0028 
6 0.0114 
fi 0.0014 
Correlation Structure Lag 
Function une tien 
16000 0.0 On0 
1 70785 Oe2 17 3 5.0000 
1.3019 0.9860 15. 0Gme 
Peole9 1.1896 20.0000 
0.2793 Hoot 25.0000 
9.4721 32.6599 25.0006 
0.2274 5.0118 30.0600 
0.0791 6.1497 45.0000 
0. 1566 4.2863 55, 0000 
0.4158 ORS Si 80.0000 
0.4850 0.9095 85.0000 
0.6269 1.4669 100.0000 
3.1748 PAB) Oi SHS. JESS) 6 OO OL0, 
176307 25.8445 150.0000 
ARS Si 25.6640 155.0000 
2.9846 20m 207 is O07 O00 
1.5724 22.036 >. OU0U 
132 22.5880 180.0000 
Ones oc 30 3a D 180.0000 
-1.6798 31.3264 20S n0000 
0.5584 27.1199 2 35 10000 
Lakes 6 22.0083 260.0000 
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TABLE C5. Mid-Intensity, Vertical Spatial Correlation 


and Structure Function 


TeastoruNc:. Variance (°C°) 


1 0.3603 x 1072 
2 0.3726 x 1072 
3 0.8536 x 10 
4 Ceaks29 
5 0.1011 
6 0.8996 x 1072 
7 0.3400 x 1072 
Correlation Structure Lag 
Pumetdon Function 
1.0000 © 0.0 ©.0 
ioe 0.2462 5.0000 
0.9265 0.3068 15.0000 
0.8702 0.1042 20.0000 
0.8345 Speen? 25.0000 
0.0002 0.0034 25.0000 
0.8742 0.4570 30.0000 
0.6962 0.3425 45.0000 
0.2549 0.7168 55.0000 
2s 0.9103 80.0000 
0.2374 1.1769 85.0000 
0.1709 0.8393 100.0000 
0.0130 0.8224 135.0000 
0.0205 Ager 7.9 150.0000 
0.0134 0.9769 155.0000 
-0.0008 0. 8464 160.0000 
61, OOS 1.3162 175.0000 
-0.0008 1.0017 180.0000 
0.0149 0.8641 180.0000 
~0.0004 0.8911 205.0000 
0.0154 0.3093 235.0000 
0.0005 0. 3356 260.0000 
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High-Intensity, Vertical Spatial Correlation 


TABLE C6. 
and Structure Function 


Thermistor No. Variance (90%) 
1 0.1800 x 107> 
2 0.2462 x 10°32 
3 0.4578 x 1072 
mn mee V0. 
5 Oneal x 10m 
6 0.5298 x 10 
i 0.2655 x 1071 
Correlation SoEuNeLvure Lag 
Function Function 
1.0000 ae (0 0.0 
0.6971 OP oSies praelelele 
0.3813 9.6565 15.0000 
0.5125 0.3485 20.0000 
1.7985 1.7245 25.0000 
Soeor 4 7 0.1730 25.0000 
We 7e2 12.0107 30.0000 
0.8521 2.9908 45.0000 
65 ES 9.2314 55.0000 
3.7301 15.2006 80.0000 
2.9532 25.7984 85.0000 
1.9051 Pome 100.0000 
0.2184 0.1375 135.0000 
0.2522 1 hoe 150.0000 
0.2551 0.6906 155.0000 
0.0156 0. 3438 160.0000 
0.0156 10.0158 175.0000 
0.0170 0.9675 180.0000 
0.3564 3.8096 180.0000 
0.0179 4.2871 205.0000 
1.0541 19.7533 235.0000 
0.0543 DAs 5) she 260.0000 
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